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ABSTRACT 

Several authors recently explored the Black-Hole mass {Mbh) vs. stellar velocity dis- 
persion (cr*) relationship for Narrow Line Seyfert 1 galaxies (NLSls). Their results are 
more or less in agreement and seem to indicate that NLSls fill the region below the fit 
obtained bv lTremaine et aT] l|20f)2f ). showing a range of cr* similar to that of Seyfert 1 
galaxies, and a lower Mbh- Until now the [O in] width has been used in place of the 
stellar velocity dispersion, but some indications begin to arise against the effective- 
ness of the gaseous k inematics in representing the bulge potential at least in NLSls. 
iBian fc Zhaol ( 2004b) stressed the urgency to produce true cr* measurements. Here we 
present new stellar velocity dispersions obtained through direct measurements of the 
Can absorption triplet (~ 8550 A) in the nuclei of 8 NLS1 gala xies. The resulting 
<7* values and a comparison with erj ml confirm our suspects (see iBotte et al.l l2004') 
that [O in] typically overestimates the stellar velocity dispersion and demonstrate that 
NLSls follow the Mbh — &* relation as Seyfert 1, quasars and non-active galaxies. 

Key words: Galaxies: active - galaxies: Seyfert - galaxies: nuclei - galaxies: kine- 
matics and dynamics 



1 INTRODUCTION 

Recently great interest was addressed to the study of 
the correlation between the mass of nuclear supermas- 
sive Black-Holes {Mbh) and the stellar velocity disper- 
sion (a*) of their hosting bulges. This relation has been 
obtain ed for nearby non-acti v e galaxies by several au- 
thors iGebhardt et a! l2000al: IFerrarese fc Merritd l2000t 
iTremaine et alJ 120021: iMarconi fc Huntll200St . supporting 
the theory that the grow th of the BH is bound to the galaxy 
formation. According to iDi Matteo et "all i2003l) this corre- 
lation, not set in primordial structures but fully established 
only at low redshift, can be understood if there is a simple 
linear relation between the total gas mass in galaxies and 
their Mbh- 

Gebhardt et al ] fcOOObl) included in their sample seven Ac- 
tive Galactic Nuclei (AGNs) and found that also these ob- 
jects followed a similar Mbh — o~* correlation of non-active 
galaxies. Other authors have then confirmed the validity of 
this result for samples of Seyfert 1 galaxies (Sis) and quasars 
<NelsorJl200d: IFerrarese et al.ll200ll: IWandeJliool iBorosonl 
l2003tlShields et all2003|)~ On the contrary, no consensus has 
been reached until now about the case of the Narrow-Line 
Seyfert 1 galaxies (NLSls). These AGNs have spectroscopic 
properties slightly different from those of the classical Sis, 
as for example, narrower permitted emission lines in the 

* Based on observations collected at Asiago observatory. 



optical wavelength domain and steeper power-law X-rays 
continua, which suggest that NLSls are hosting BHs with 
smaller masses accreting at high rates, close to the Edding- 
ton limit. This hypotesis would imply that NLSls have also 
less massive bulges than classical Sis, of course assuming 
that their Mbh are correlate d to the physical p r operties 
of the ir ho sting bulges. In deed, IWang fc Lul feOOllhlWandel 
J2002I) and lWandel ( 2004) found that there is no clear differ- 
ence in the Mbh — cr* relation between narrow-line, broad- 
line AGNs and non-acti ve nearby galaxies, while opposite 
result s were obtained bvlMath ur 1 Kuraszkiewiczj fc C zernvl 
j200ll).lBian fc Zhaol i2004al) . iGrupe fc Mathurl <2004l) . and 
IBotte et all J2004I) . The main problem is that cr* in AGNs 
is available f or few Sis and only 3 N LSls: Mrk 110 (90 ± 7 
km s" 1 ) b y IFerrarese et alJ < |200lt> . NGC 4051 (88 ± 13 
km s _1 ) bv lNelson fc Whittle! ]l995l). an d Mrk 766 (106L40 
km s _1 ) bv |jmil§n^BernTo^taL| J2000I) . In place of <r», the 
[O ill] emission linewidth was extensively used as a represen- 
tation of the bulge velocity dispersion. 

Here we present new direct determinations of <r» in the 
nuclei of 8 NLS1 galaxies, obtained using the Call absorp- 
tion triplet (8498. 8542 8662 A). The targets were extracted 
from IBotte et all (I2004T) . except for Mrk 766. In Section 2 
we present the spectroscopic observations and the data re- 
duction; in Section 3 we give a detailed description of the 
methods we used to measure stellar and gaseous kinemat- 
ics, and we present the results. Finally, in Section 4 we dis- 
cuss the implications of our results on the validity of the 
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Figure 1. The rest-frame spectra of the 8 NLS1 galaxies showing 
11559) are also plotted for comparison. 

FWHM([Om])/2.35=cr» and M BH - a, relations for NLS1 
galaxies. 



2 OBSERVATIONS AND DATA REDUCTION 

We have observed the nuclei of 7 NLSls with the Asiago 
Faint Object Spectrograph Camera (AFOSC) mounted at 
the 1.82 m telescope of the Padova Astronomical Obser- 
vatory (Asiago, Italy). A Volume Phase Holographic grism 
was used in combination with a 1.69 arcsec-slit to cover the 
region ~ 8200-9200 A which includes the Call triplet ab- 
sorption lines, with a dispersion of ~ 0.88 A/px and a rela- 
tively high spectral resolution (~ 50 km s^ 1 ). Mrk 766 was 
also observed in the optical wavelength range (4300-8000 
A) at medium resolution, in order to d erive the Mbh and 
[O ill] width information not availble in lBotte et alJ (2004) 
and useful for the following analysis. In addition to galaxy 
observations, the spectra of two template stars, HR 294 (G9 
III) and HD 11559 (K0 III), were obtained for kinematic 
analysis with the same configuration. Two spectral types 
were chosen to minimize the template mismatch effects. 

Mrk 896 was extracted from the Isaac Newton Group 
(INC) Archive. A spectrum covering a larger range (~ 7700- 
9400 A) was obtained in August 1994 at the William Her- 
schel Telescope (Canary Islands, Spain) with ISIS Double 
Beam Spectrograph. The R316R grating used in combina- 
tion with a 1.47 arcsec-slit yielded a dispersion of ~ 1.35 
A/px and an instrumental resolution of ~ 65 km s~ . A 
standard star (BD+17 4708) was observed during the same 
run for spectrophotometric calibration. 

The spectra are shown in Fig.Q Details about the ob- 
servations are given in Tabled 

All data were reduced with the same procedure. The usual 
reduction steps, that is bias and flat field corrections, cosmic 
rays removal, wavelength calibration by means of compari- 
son lamps and sky-background subtraction, were carried out 
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e Ca II triplet absorption lines. The template stars (HR 294 and HD 



Table 1. Observation Log. 



Galaxy z AA Date Texp 

(A) (sec) 



Mrk 359 





017 


8200- 


-9200 


2003- 


09 


08 


7200 


Mrk 335 





026 


8200- 


-9200 


2003- 


12 


07 


9000 


UGC 3478 





013 


8200- 


-9200 


2003- 


12 


08 


7200 


Mrk 705 





029 


8200- 


-9200 


2003- 


12 


08 


7200 


IR 04312 





020 


8200- 


-9200 


2004- 


01 


26 


7200 


IC 3599 





022 


8200- 


-9200 


2004- 


01 


26 


7200 


Mrk 766 
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-6500 


2004- 
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2400 








6200- 


-8000 
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1800 
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7700- 


-9500 
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08 


16 


1800 



with IRAF packages 1 . The optical spectrum of Mrk 766 was 
also flux calibrated through the observation of two spec- 
trophotometric standard stars, FEIGE56 and G191-B2B. 

For each object we extracted the nuclear spectrum sum- 
ming a number of pixels along the slit to increase the signal- 
to- noise ratio (S/N). An aperture of ~ 3.3 arcsec (7 px) was 
chosen for the Asiago data on the basis of the seeing, which 
ranged between 2 and 3 arcsec during the observations. This 
aperture corresponds to ~ 0.6 — 2 kpc for galaxies in the 
range 0.01 < z < 0.03. For the ING spectrum we used an 
aperture of 1.7 arcsec (5 px), corresponding to ~ 0.9 kpc, 
following similar arguments. 



1 IRAF is written and supported by NOAO (Tucson, Arizona), 
which is operated by AURA, Inc. under cooperative agreement 
with the National Science Foundation 



3 DATA ANALYSIS 

3.1 Stellar velocity dispersion 

To measure a* we have used the Fourier cross-correlation 
method as implemented i n the IRAF task FXC OR which is 
based on the method of iTonrv fc Davis! il979t) . The spec- 
trum of the galaxy G and the stellar template T are resam- 
pled into N bins, where each bin number n is proportional 
to In A, and cross-correlated in Fourier space: 

/oo 
T(k)G{k + n)dk (1) 
-oo 

where x means cross correlation. 

This method assumes that the G(n) is a convolution of 
T(n), representing the true stellar population in the galaxy, 
with a broadening function b(n), offset from zero velocity 
by an amount 5: 

/oo 
b(x — S)t(n — x)dx (2) 
- oo 

where * denotes convolution. 

The maximum peak of the cross-correlation function 
(CCF) is fitted by a smooth symmetric function (in general 
a gaussian function) , whose location and width are related to 
the galaxy redshift and velocity dispersion of the stars in the 
galaxy. In order to recover a* from the F WHM of the CCF 
peak, we adopted the suggestions given bv lNelson fc Whittle! 
(1995). We convolved the template spectrum with gaus- 
sian functions of known increasing dispersions (from 10 to 
300 km s" 1 ) and we measured the corresponding increas- 
ing widths of the CCF peak. Then, we fitted the empirical 
FWHM-tr, relationship with a polynomial function. 

Before to run FXCOR, the spectra of the NLSls were 
moved to the rest-frame and their continuum carefully re- 
moved. This last step proved to be rather important, since 
the cross-correlation method is sensible to low-frequency 
fluctuations and our targets are often characterized by emis- 
sion lines and/or bumps. In particular, 5 out of the 8 ob- 
served targets show the O I A8446 emission line, and in at 
least 3 galaxies (Mrk 335, Mrk 896, UGC 3478) the Can ab- 
sorptions are overlapped to broader Can emissions. Then, 
FXCOR was run using two template stars for the Asiago 
data, and the only one available for the ING archival data. 
The resulti ng CCF peaks had va lues of the Tonry-Davis R 
parameter iTonrv fc Davisl|l979l) . a sort of S/N indicator, 
ranging from 8 to 20. No filtering of high-frequency noise 
was applied, since we verified that it was not necessary and 
introduced variations of the results depending on the arbi- 
trary choice of the filter shape. 

To esti mate the uncer t ainty of the cr* measurements, 
we followed lOhvama et al.l (I2002T) . Firstly, two spectral re- 
gions were independently correlated for each template spec- 
trum: 8300-8800 A and 8450-8750 A. This allowed to 
have four values of <r* for each galaxy. We adopted the 
mean of these measurements and calculated the RMS (Sai). 
Sai was not estimated for Mrk 896, since only one tem- 
plate star was available. Then, we made 20 spectra for 
each object by adding artificial noise with the IRAF task 
MKNOISE to increase the RMS by ~ 20 per cent. FX- 
COR was run again 20 times with these new spectra and 
the scatter of the peak width measurements was adopted 
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Table 2. Stellar and gaseous velocity dispersions, and BH masses. 



Galaxy 




CT [0 III] 






(km s _1 ) 


(km s _1 ) 


fxlO 6 Mr-,) 


Mrk 335 


64 ± 17 


140.8 ± 3.4 


8.61 


Mrk 359 


112 ± 12 


67.0 ± 0.4 


0.69 


UGC 3478 


89 ± 18 


109.8 ± 2.6 


0.81 


Mrk 705 


82 ± 25 


185.9 ± 2.1 


46.13 


IR0 4312 


64 ± 20 


136.8 ± 7.8 


0.83 


IC 3599 


85 ± 17 


106.9 ± 4.4 


0.13 


Mrk 766 


81 ± 17 


131.2 ± 0.7 


0.63 


Mrk 896 


87 ± 11 


129.6 ± 0.6 


4.43 


Mrk 110 


90 ± 7 a 


123.3 6 


6.65 c 


NGC 4051 


88 ± 13 d 


81 b 


1.35 c 



a Average value from Fcrrarese^^^L (2001) 
b From Whittle ( 19921 
c Average value from Kaspi et al. (2000) 
d From Nelson fc Whittle (1995) 



as error of the CCF peak fitting (<5<T2). Finally, the un- 
certainty of the instrumental resolution (<5r ~2-4 km s _1 ) 
was consider ed as additional s ource of error by deriving the 
Eq. (18) bv ITonrv fc Davis! Jl979t) : Sa 3 = 2t8t/g,. As- 
suming these errors are independent one from each other, 
the total Atr, was obtained throug h the following relation: 
Act, = yJ{8(Ji) 2 + (8a 2 ) 2 + (Sa 3 ) 2 . The stellar velocity dis- 
persions are listed in Tabled 

3.2 Gas velocity dispersion 

In addition to measurements of stellar kinemat ics from ab- 
sorpt ion lines, we took from our previous work jBotte et alJ 
2004) the optical spectra of the same objects, containing 
the Narrow Line Region (NLR) emission lines which allow 
to measure the gas kinematics. In particular, the width of 
[O ill] A5007 line was measured and then corrected for in- 
strumental line broadening to determine the velocity dis- 
persion of the gas in the NLR. Since the [O ill] emission line 
is bright in these spectra, its profile is very little affected 
by noise. Therefore, we fitted [O ill] with gaussian functions 
only changing for 5 times the continuum level, chosen by vi- 
sual inspection. The average dispersion <T[o hi] and its rela- 
tive scatter were calculated for each galaxy. The gas velocity 
dispersions are given in Table [5] 

3.3 Results 

In Fig.|5|we plotted cr* values vs. [O ill] widths (<T[o hi]). The 
result is very interesting since it clearly shows that [O ill] 
sistematically overestimate the stellar velocity dispersion in 
NLSls. In fact all targets except one (Mrk 359) occupy the 
upper half of the a, — ayo in] plane, and are not grouped 
around the 1:1 separation line. I n other words, this is the 
first indication that using ayo nil in place of a* could be 
wrong in NLSls, which seem to have a < a, >~ 84 ± 14 
km s~ . In this picture, Mrk 359 is an outlier. It is the only 
case among o ur data having 070 rr n < A justification 
was given bv Ijimenez-Benito et al.l i2000Tl . who suggested 
that such cases could be caused by gas settled in a cold 
rotating disc. 
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Figure 2. Stellar vs. gas kinematics. Solid triangles are the ob- 
served galaxies, while open triangles the two additional NLSls 
taken from literature. The solid line is o^q jjji = cr, . 



Figure 3. Stellar velocity dispersion vs. BH mass. The solid line 
is the lTremaine et alj i2002fl fit. Symbols are like in Fig. |3] 



Fig.|2]shows the new a, data plotte d vs. the BH masses . 
Mbh values (Table are taken from iBotte et al] (l2004f) 
except for Mrk 766, which was observed also in the optical 
range between 4000 and 6000 A to measure the FWHM(H/3) 
and the continuum luminosity at rest-frame A = 5100 A 
as well as the [O ill] width. Its <j[o m] and BH mass value 
are also listed in Table |2] In both figures, open triangles 
indicate Mrk 110 and NGC 4051, two additional NLSls 
with known a,, mentioned in Section 1. Their <T[o hi] and 
Mbh values are taken from literature. For Mrk 766 we used 
our measurements, which a re however in agreement with 
Ijimenez-Benito et al] i2000ll within the error bars. 

It is clear from Fig. [3] that th ese new a, values fill a 
range more in agreement with the iTremaine et al] (I2002T) 
relation: log Mbh = 4.02 log cr* — 1.12, considering a mass 
range of ~ 10 6 - 10 7 Mq for BHs in NLS1 galaxies. In fact, 
after having plotted our data over the Mbh — o* relation, we 
see that NLSls are well arranged around the Tremaine et al. 
fit, and not sistematically below the fit as previously found. 
The data show a significant scatter around the fit, not very 
different from that observed in ordinary Seyfert 1 galaxie s 
(see e.g. lWang fc Lj200ltlWandei200alBotte et aljEool . 
Moreover, 5 out of 10 NLSls have a factor of 3-10 lower BH 
masses than predicted from the Mgii — <r* relation. This 
would suggest that NLSls have a tendency to lower Msg/a, 
ratios relative to non-active galaxies and broad-line AGNs, 
but the limited number of points does not allow to draw a 
conclusion about this point. 



4 DISCUSSION AND CONCLUSIONS 

Several authors tested the Mbh — o"* correlation validity 
for AGNs assuming that the gaseous component could rea- 
sonably trace the nuclear gravitational potential. In partic- 
ular, the [O in] width was used as surrogate of the stellar 
velocity dispersion, essentially for a practical reason: the 
[O ill] emission line is usually bright in AGNs and there- 
fore easily measurable even in case of distant objects. On 
the contrary stellar absorptions require longer integration 
time and/or large telescopes to reach the S/N ratio neces- 
sary to give a reliable estimate of the stellar velocity dis- 
pe rsion. The FW HM([Q ill]) to u» conversion was proposed 
bv lNelson fc Whittle! ]l996l) . who observed a large sample 
of AGNs in two spectral ranges, one including H/3, [O ill] 
emission lines and Mg b triplet absorption, and the other 
containing the Call triplet absorption. The data of 66 ob- 
jects were plotted showing a m oderate correlation with a sig- 
nificant scatter. More recently I Jimenez-Benito et al] <|2000T) 
obtained similar results on a sample of Seyfert 1 and 2 galax- 
ies, demonstrating that the majority of the galaxies form a 
cloud of points around the 1:1 line of the gas vs. stellar ve- 
locity dispersion plot, with very weak or even no correlation. 

None the less, the use of the [O in] width i n plac e 
of direct a, measure ments led iBian fc Zhaol ]2004aT) , 
iGruoe fc Mathurl ]2004l) . and lBotte et all J2004ir to demon- 
strate that Sey fert 1 galaxies agree w ith the Mbh — o* re- 
lation given bv ITremaine et al J J2002T) . even if with a large 
scatter of cr, (from 50 to more than 500 km s _1 ), while 
NLS1 galaxies do not. Moreover, NLSls appear sistemati- 
cally arranged below the Tremaine et al. fit. Interestingly, 
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this result was already included in lBorosonl i2003f) . who car- 
ried out an analysis of the Mbh~[0 hi] A5007 width relation 
for 107 low-redshift radio-quiet QSOs and Seyfert 1 galax- 
ies and found that this correlation is real, but the scatter is 
large, about 0.2 dex. Really, it is easy to verify that the nar- 
row line AGNs (FWHM(H/3) < 2000 km s _1 ) of his sample 
are mostly located below the Tremaine et al. fit. Therefore, 
authors were induced to believe that NLSls do not follow 
the BH-bulge correlation as Seye fert 1, quasars an d normal 
galaxies. But, very recently. iBian fe Zhaol J2004bl) obtained 
a good consistency in BH mass estimates of NLS1 galax- 
ies by using H/9 and soft X-ray luminosity, while no agree- 
ment was found with the [O ill] linewidth. Their result casts 
doubt on the reliability of [O ill] as indicator of stellar ve- 
locity dispersion, at least for NLSls. And, indeed it should 
be mention ed that the FWHMffO m1 )/2.35=cr« relation was 
obtained bv lNelson fe Whittle! <ll996t) for a sample of Seyfert 
1 and 2 galaxies, but it was never verified for NLSls. 

The situation clearly required new direct cr» measure- 
ments. We obtained these values for 8 NLSls observing their 
Can triplet. The choice of this spectral range was induced 
by the fact that NLSls are very often characterized by large 
Fen emission multiplets (4400-4700 A and 5150-5350 A) in 
correspondence of the absorption lines commonly used to es- 
timate a* , like Mg I A517 5 and Fe I A5269. Unfortunately, as 
shown bv lPerssonl il988f) . NLSls with strong Fell are likely 
to have Ca II in emission rather than in absorption. 

By comparing stellar and gaseous kinematics we found 
that the FWHM([Om]) /2.35=<7„ relation does not seem 
to hold for NLSls, which are arranged mostly at values of 
°"[o in] larger than expected. Contrary to the recent find- 
ings, the direct a, measurements span a narrower range 
typically lower than 100 km s _1 , and such range is in agree- 
ment with the stellar veloci ty dispersions expected by the 
Mbh — cr* relation given bv lTremaine et al] B2002I) in case 
of M BH ~ 10 6 - 10 7 M Q . Indeed, with our cr* determina- 
tions, NLSls are arranged around this relationship, suggest- 
ing that in NLS1 galaxies BH and bulge correlate, and lower 
mass BHs correspond to lower stellar velocity dispersions, or 
in other words lower mass bulges. We st ress that this resul t 
is in agreement with our previous work jBotte et aljl2004) . 
where we showed that NLSls are mostly confined to the 
lower ranges of the MBH-'Lbuige(B) plane, and are well cor- 
related to Seyfert 1 galaxies and quasars. The fact that the 
stellar velocity dispersion in place of the gaseous one yields a 
tighter relation when correlated with BH mass, is essentially 
caused by the star kinematics which is more representative 
of the bulge gravitational potential than the gas in AGNs. 
Even if the gas motion is mainly controlled by the mass of 
the bulge, several observations in the past showed that the 
base of the line profiles emitted by the Narrow Line Region 
(NLR) in many active galaxies is hardly reproduced sim- 
ply assumin g a pure virial motion in the potential of the 
host galaxy iVeilleuxlll99ll . and references therein) . Other 
important factors or processes can be involved, and these 
additional processes can be at the origin of the larger scat- 
ter observed when the FWHM([0 III]) is used as surrogate 
of a,. The interaction bewteen the gas and the ejected ra- 
dio plasma observed in many Seyfert galaxies may strongly 
influence the gas kinematics and contribute with o utflow- 
ing motions of the NLR clouds jNelson fe Whittle! Il99r3) . 
Moreover, supersonic winds generated by the active nucleus 



are believed to accelerate the smaller clouds of the NLR lo- 
cated closer to the central engin e and sometimes dominate 
their virial velocity dSmithlll993D . This effect could be more 
important in NLS1 galaxies, since stron g nuclear winds are 
expec ted in case of high accretion rates (|Leighkv_& Halpern 
2000). In addition, it was pointed out bv TNelson^^VhrtHe 
(1996) that a slight tendency exists for barred and/or dis- 
turbed active galaxies to have broader [O III] emission lines. 
In conclusion, we obtained results with important implica- 
tions, which strongly encourage additional observations to 
enlarge the sample of the observed galaxies. 
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